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The analysis of the EM radiation from a single charge shows that the radiated 
power depends on the retarded acceleration of the charge. Therefore for consis- 
tency, an accelerated charge, free from the influence of external forces, should 
gradually lose its acceleration, until its total energy is radiated. Calculations 
show that the self force of a charge, which compensates for its radiation, is 
proportional to the derivative of the acceleration. However, when using this 
self force in the equation of motion of the charge, one gets a diverging solution, 
for which the acceleration runs away to infinity. This means that there is an 
inconsistency in the solution of the single charge problem. However, in the con- 
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struction of the conserved Maxwell charge density, there is implicitly an integral 
over the corresponding world line which corresponds to a collection of charged 
spacetime events. One may therefore consistently think of the "self force" as 
the force on a charge due to another charge at the retarded position. From 
this point of view, the energy is evidently conserved and the radiation process 
appears as an absorbing resistance to the feeding source. The purpose of this 
work is to learn about the behavior of single charges from the behavior of a 
real current, corresponding to the set of charges moving on a world line, and 
to study the analog of the self force of a charge associated with the radiation 
resistance of a continuum of charges. 
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1. Introduction 

Antennas and EM radiation problems have been extensively discussed in the 
literature, and usually have closed form solutions, or at least closed numerically 
formulated solutions [1,2]. This means that the behavior of charges which are 
part of a real given electrical current is well described. 

However, the behavior of single charges has still many open questions [3,4,5]. 
In the construction of the conserved Maxwell charge density, there is implicitly 
an integral [6] over the corresponding world line which corresponds to a col- 
lection of charged spacetime events. One may therefore consistently think of 
the "self force" as the force on a charge due to another charge at the retarded 
position. In this paper we wish to deduce and extrapolate the behavior of single 
charges, from the behavior of a real electrical current. Specifically, we would like 
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to understand if the " self force of a charge, which acts as a radiation damping 
force, is an analogy to the radiation resistance, which acts like a power absorber 
from the feeding source. 

We will use in this paper e = l/(47r);^ = An; c = 1, and put back the 
standard values only where it will be specifically mentioned. 

2. Formulation 

We formulate here the following problem: 

z 

L/2 



-L/2 

Figure 1: Short dipole antenna 

The antenna is fed by the current source I(t) (at z = 0), which can be ex- 
panded in a Fourier series of functions. Therefore, without loss of generality we 
may consider the harmonic excitation: I(t) = I cos(u;t). At an arbitrary z, the 
current will be: 

I{z,t) = I f (z,t) + I b (z,t) (1) 
where //, h are the current waves which move forwards and backwards, respec- 



io cos(wt) 
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tively, in the z direction, and have the form of: 



If,b{z,t) = Ci,2COSw(t=FZ + 01,2) 



(2) 



where the constants Ci t 2, <t>i,2 must satisfy the boundary conditions of zero cur- 
rent at z = ±L/2 and Iocos(ujt) current at z = 0. Those boundary conditions 
establish the spatial dependency of the current f(z) so that the current in (1) 
may also be written as: 



However, these boundary conditions are not important for the work of this pa- 
per because we are interested in locally analyzing the behavior of charges. But 
we will use the general form of the current, for establishing connections between 
the current and microscopic parameters, like speed and acceleration of the single 
charges. So we rewrite the current in (1) as: 

I(z, t) — C\ cosuj(t — z + <j>i) + C 2 cosuj(t + z + 4>2) (3) 



where / is the frequency of the excitation and to is the wavelength. 

Next we have to calculate the speed of the charges v(z,t). Having a one 
dimensional configuration, it will be convenient to define p, the charge density, 
as the charge per unit of length. Therefore, to be consistent we have to define 



I(z,t) = I cos(ut)f(z), /(0) 



= l,/(±L/2) = 



(la) 



and 



lo = 2nf = 2n/X 



(4) 
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J, the current density, as the current I itself: 

J(z, t) = If(z, t) + Ib(z, t) = C\ cosuj(t — z + 4>i) + C2 cos u)(t + z + fo) (5) 
Knowing that: 

3^ = therefore + |f = (6) 
we calculate: 

|f = w[-Ci sinw(t - z + <j>i) + C 2 sinw(t + z + <j> 2 )] (7) 
Therefore: 

p = uj J dt[— C\ smuj(t — z + 4>i) + C 2 sino;(t + z + (f> 2 )] = 
I f (z,t)-h(z,t)+p (8) 

po is the integration constant and represents the density of the free electrons 
per unit of length in the conductor. As we shall see, p S> \If(z,t) — I b (z,t)\ 
therefore If(z,t), Ib(z,t) constitute low amplitude waves moving forward and 
backwards, respectively, around the average charge density po- 
Knowing the relation J = pv, we obtain for v(z, t): 

V ( Z +) = 1 = L = If(zS)+h(zt) , s 

u \ z i p p i f (z,t)-h(z,t)+p y y > 

For any practical purpose po is the dominant term in the denominator of 
(9), since p 3> Ci^/c (we insert here the light velocity c). For example, if 
Ci.2 — 10A, then Ci^/c = 1.87 x 10 11 electrons/m. Let's say our wire has a 
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surface of 0.5mm 2 , so Ci^/c corresponds to 7.5 x 10 17 electrons/m 3 , while the 
charge density in copper is 8 x 10 28 electrons/m 3 . So p is bigger than Ci^/c 
by 11 orders of magnitude, and the velocity of the electrons is smaller than c by 
11 orders of magnitude, i.e., about 10~ 4 m/sec. We therefore may approximate 
the velocity of the charges by: 



v(z,t) = j- o I(z,t) 



(10) 



3. Discretization of the continuum 

The next step is to discretize the current into charges of magnitude, having a 
distanceone from the other, so that, as in Figure 2. 



conductor 



Charge A 



Charge B 



free charges 
'of size q 



~f 

Az = q/p 
_t 



Figure 2: Close-up view of the conductor 

We will examine the forces acting on the "test charge" B (Fig. 2), as a result 
of a disturbance produced on charge A. 

Before the disturbance occurred, all charges had a constant velocity, and on 
each charge, two opposite static forces of magnitude q 2 /(Az) 2 acted. 
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A disturbance which occurred at charge A at an earlier time t = — At = 
—Az, appears as a disturbance perceived at charge B at t = 0. 

The motion of the "disturbed" charge A, z A (t) can be expanded for t = —At, 
up to third order, as follows: 

z A (-At) = z Ao - vAt + \aAt 2 - ±dAi 3 (11) 

where v, a and a are the velocity, acceleration and its derivative, measured at 
t = 0. The velocity of the "disturbed" charge A is: 



v A (-At) = v - aAt + \aAt 2 (12) 

The full distance to zb from the disturbed charge at A must arrive at light 
velocity, therefore, Az at time — At is: 



Az(-At) = z B - z A (-At) = Az a + vAt - ±aAt 2 + idAt 3 = At (13) 

Formula (13) tells us that the retarded distance equals the time needed to prop- 
agate the disturbance, and expands this distance in powers of At, where Az is 
the distance between the charges A and B at t = 0. 
The near field of a charge is expressed by [4]: 



E * ~ IS^Fret - [ 7 A*(i-„ A )P|_ at > ( 14 ) 



where T is very close to 1. We will use (12) and (13) for putting the values of 
Az and v A into (14) and get from (13): 
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Az - At(l -v+ \aAt - iaAt 2 ) = 



(15) 



With the aid of (15), we obtain from (12) 

l-v A {-At) = ^ -±aAi+±dAt 2 + aAi-±dAi 2 ~ l + ±aAt-±dAi 2 (16) 

Here, we have used the approximation Azo — At because the velocity is very 
small. 

Using this result, Az = At and 7 ~ 1 in (14), we get the field at charge B, 
caused by a disturbance which occurred At earlier at charge A, as 



The first term q/At 2 is always cancelled by the force of the "other" neighbor, 
as mentioned before, it therefore can be ignored. 

The third term is identical to what is considered to be the field which creates 
the self-force of a charge, but here it was derived as the force on a charge, due 
to a disturbance on a neighboring charge. As we shall see, it is the only term 
which creates the radiation resistance (which we see is a local phenomenon on 
the world line). 

We will call the last 3 terms of (17) E^amp, and after replacing At by Az, 
we obtain: 



The potential difference V on a wire segment of length Az, resulting from E^amp, 
will be denoted as the "damping" tension (or "damping" voltage) and is calcu- 



(17) 



Edamp — <l(—-ztz + 3^ ~ I° 2 ) 



(18) 
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lated by: 

V = -E damp Az (19) 
The magnitude Az is according to the discretization defined above, so that 

Az = q/p (20) 
We therefore get from (18), (19) and (20): 

v = _ q 2 [2/3v-i,mf] +qi (21) 

According to (10): 

v = ±§ and v = ±%l (22) 

po at pa at 2 \ > 

Dealing with harmonic excitation, d/dt is like multiplication by to (up to a 90° 
phase). So v ~ ui 2 v, and (v) 2 ~ (wji) 2 . Having u smaller by 11 orders of 
magnitude than light velocity, v 2 is completely negligible relative to v. 
We therefore may write (21) (after neglecting the v 2 term), using (22): 

y = -Uf a ) 29 ^ + = -\^ 2 U + (23) 

The power radiated by the segment Az is AP = VI and given by: 

AP=-I^A^ + f/ A , (24) 
According to (3), (and also because I satisfies the wave equation): 
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d 2 I _ , ,2r _ /'ZTr^r 



(25) 



The constant ratio between d 2 I/dt 2 and I means that the current is "in phase" 
with its second derivative, and therefore the first part of AP in (24), integrated 
over time, represents radiated energy. However, the multiplication of I by dl jdt 
has the form of cosw(t — z + 4>) sin w(i — z + 4>) and therefore the second part of 
AP represents a reactive power, which results in zero energy after integrating on 
an integer number of time cycles. The reactive power represents power which is 
returned to the source each time cycle. We are therefore interested in the first 
term of AP in (24). 

Putting (25) into the first part of (24) we obtain: 



AP(z) = f(f) 2 I[/(z)Az] 2 (26) 



In (26) we omitted the time dependence, which is harmonic (see (la)), and 
multiplied by \ which is (1/T) ^ cos 1 {2-k / T)tdt . So (26) represents the mean 
emitted power over a time cycle T = 2-k/uli. 

Therefore the mean power emitted from the whole wire is: 

P= im 2 W%I(z,t)dz? ee |(f ) 2 |[7 „L] 2 = | ( 2.) 2 (^) 2 i4 2 , (27) 

Here we defined I av ee (1/L) /ff /2 I(z)dz. 
The radiation resistance is defined by: 

P^ = i(2^m) 2 (W/ ) 2 . (28) 
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The factor 2 is due to P being averaged on a time cycle, and Iq is the ampli- 



impcdance of the free space, if we had used eo = 1; but we used e = 1/4-7T, so 
we have to multiply (28) by Z = = 30. We therefore obtain the radiation 
resistance: 



which is identical to the result in the literature [1,2]. 
4. Discussion 

The radiation resistance of two particular cases can be immediately derived 
from (29). 

One is a halfwave dipole L = A/2, for which the boundary conditions give 
f(z) — cos(27rz/A) in (la). Here we obtain I av = (2/tt)Iq, hence R — 80f2. 

The second is an ideal case of very short dipole antenna, in which we suppose 
lav = h , (as if the boundary conditions were satisfied by two discs at the ends of 
the wire to act as charge reservoirs). In this case we obtain R — 807r 2 (L/A) 2 fl 
This result could have been obtained by directly dividing the voltage in the first 
part of (23) by the current I = I , and obtain R = 80tt 2 (Az/A) 2 ^, if L = Az. 

This is actually the most interesting case for our study, because it considers 
the segment Az (which contains the charge q) as an antenna. This segment has 
the radiation resistance R = (2/3)(27r) 2 (Az/A) 2 (in our units), which directly 
results from the third term of (17), and could also be written as: 



tude. Now we put back the units: Z a = l/(e c) = ^/^ / e o = 1207T would be the 



R = 80n 2 (j;) 2 (I av /I ) 2 n 



(29) 



R = 




(30) 
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where the work is P se if = F se ifV, F se if being (2/3)q 2 a, and I is the charge flow 
per unit time - or current. We see that F se if in a continuum of charges could 
be considered as the force acted by the perturbed neighbor. 

The radiation resistance is therefore the "self work" divided by the charges 
flow rate squared. 

We have to remark that this result has been derived with the aid of several 
approximations like pa 3> Ci^/c, resulting in v <C c, which are very accurate for 
charges in a real wire, but cannot be counted on for a single charge. Therefore, 
result (30) is not a proof of the formula F se if — (2/3)q 2 a, but rather shows that 
in a macro perspective, it looks like the self force of a charge, which accounts 
for the radiation, is (2/3)q 2 a. 

Moreover, P se if comes out to be also (2/3)q 2 a 2 (averaged on a time cy- 
cle), which is the Larmor formula, describing the emitted energy per unit 
of emission time. This is because: v(t) — = 7 " cos ^ t ) = Vq cos(wi), 

w Po Pa u v /j 

a(t) = — covo sin(wt), a(t) = — lo 2 vq cos(wt) and therefore dv = — (luvq) 2 cos 2 (tut) 
and a 2 = (ujvo) 2 sin 2 (wt). Both sin 2 (o;t) and cos 2 (wt) averaged on a cycle, give 
1/2. 

The result that \av\ = a 2 implies that P 8e if represents the "self work", i.e. 
F se if v and also the emitted energy per unit of emission time, suggesting that 
the self force is an energy conservation force. 
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